Energy  Conversion  and  Management  86  (2014)  992-1000 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Energy  Conversion  and  Management 

journal  homepage:  www.elsevier.com/locate/enconman 


Energy 


Conversion 

IManagement 


Simulation  and  evaluation  of  a  CCHP  system  with  exhaust  gas 
deep-recovery  and  thermoelectric  generator 

J.L.  Wang,  J.Y.  Wu*.  C.Y.  Zheng 

Institute  of  Refrigeration  and  Cryogenics,  Shanghai  Jiao  Tong  University,  Dongchuan  Road  800,  Shanghai  200240,  China 


CrossMark 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  26  March  2014 
Accepted  12  June  2014 
Available  online  10  July  2014 


Keywords: 

CCHP 

Deep-recovery 

TEG 

Exhaust  gas 
Simulation 


Combined  cooling,  heating  and  power  (CCHP)  systems  are  thought  to  be  highly  efficient  in  energy 
utilization.  But  there  are  still  potentials  to  further  improve  system  performance.  This  work  proposed  a 
CCHP  system  based  on  internal  combustion  engine  (ICE)  for  power  generation,  refrigeration  and  domestic 
hot  water  production.  Thermoelectric  generator  (TEG)  and  condensing  heat  exchanger  are  applied  to 
efficiently  recover  the  exhaust  gas  waste  heat  of  ICE.  All  the  energy  flows  are  designed  based  on  energy 
cascading  utilization  principle.  Basing  on  the  test  results  of  a  16  kW  ICE,  CCHP  system  characteristics  are 
investigated  by  simulation  from  idling  condition  to  full  load  condition.  Especially,  the  part  load  perfor¬ 
mance  of  TEG  and  absorption  chiller  are  simulated  and  discussed.  The  feasible  operating  zone  of  ICE 
and  feed  water  flow  rate  are  figured  out  to  keep  the  domestic  hot  water  temperature  within  a  certain 
range.  Results  show  that  the  primary  energy  efficiency  of  system  can  reach  0.944,  thanks  to  the  condens¬ 
ing  heat  recovery  from  exhaust  gas.  The  primary  energy  saving  ratio  and  cost  saving  ratio  can  reach  0.304 
and  0.417,  respectively.  Considering  some  more  equipment  is  incorporated,  the  total  investment 
increment  is  about  11.1%. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  application  of  combined  heating  and  power  has  a  history  of 
more  than  100  years,  and  it  has  grown  into  combined  cooling, 
heating  and  power  (CCHP)  technology  which  is  thought  as  a  most 
advantageous  branch  of  distributed  generation  [1-3].  Wu  and 
Wang  [3]  and  Caresana  et  al.  [4]  report  that  in  CCHP  system,  espe¬ 
cially  in  micro  and  small  size  system,  internal  combustion  engine 
(ICE)  is  the  most  suitable  prime  mover,  for  its  good  adaptability, 
good  reliability,  and  low  initial  investment.  The  efficiency  of  an 
ICE  is  usually  inferior  to  that  of  a  centralized  power  plant,  even  if 
considering  the  transmission  loss  of  public  grid.  But  the  utilization 
of  waste  heat  from  ICE  can  make  it  be  advantageous.  So,  it  is  quite 
significant  to  make  the  most  of  waste  heat  from  ICE.  At  present,  the 
reported  primary  energy  efficiency  of  commonly  used  CCHP  sys¬ 
tems  ranges  from  0.5  to  0.86  3,5-8].  Technically,  there  are  poten¬ 
tials  to  further  improve  the  performance  of  CCHP  system. 

Firstly,  high  efficient  thermal  activated  devices  should  be  used 
as  much  as  possible.  The  exhaust  gas  temperature  of  ICE  can  be 
over  500  °C  [5].  When  there  is  cooling  load,  the  exhaust  gas  can 
be  used  to  power  a  double  effect  absorption  chiller  (DAC)  directly. 
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In  many  CCHP  systems  [5,9-12],  exhaust  gas  is  recovered  by  hot 
water,  and  then  the  hot  water  is  used  to  power  a  single  effect 
absorption  chiller  (COP  of  about  0.7  [13-15]),  or  an  adsorption  chil¬ 
ler  (COP  of  0.4-0. 5  [13,16]).  The  COP  of  DAC  powered  by  exhaust 
gas  can  range  from  1.2  to  1.48  [15,17-19],  which  is  obviously  much 
higher  efficient  than  other  commonly  used  thermal  activated 
devices.  So  DAC  should  be  considered  as  much  as  possible. 

Secondly,  the  exhaust  gas  can  be  used  to  power  thermoelectric 
generator  (TEG)  for  electricity  generation  instead  of  just  heating 
hot  water.  The  quality  of  electricity  is  much  higher  than  that  of 
hot  water.  Organic  Rankine  cycle  is  also  a  good  choice  to  recover 
the  waste  heat  of  exhaust  gas  [20].  But  the  additional  system  is 
far  more  complex  than  TEG.  Bass  et  al.  have  developed  a  TEG, 
which  was  installed  outside  the  exhaust  gas  pipe  of  a  truck  and  a 
capability  of  1  kW  power  output  was  demonstrated  [21,22].  Ikoma 
et  al.  [23]  and  Karri  et  al.  [24]  have  also  respectively  studied  TEG 
systems  applied  in  exhaust  gas  waste  heat  recovery  of  engines.  A 
TEG  can  be  compactly  integrated  on  the  outside  of  exhaust  gas 
pipe,  and  it  could  be  feasible  in  a  CCHP  system. 

Thirdly,  exhaust  gas  deep-recovery  can  be  applied  to  make  the 
most  of  waste  heat.  The  exhaust  gas  contains  much  gaseous  water, 
especially  for  combustion  of  natural  gas  (NG).  If  the  exhaust  gas 
can  be  cooled  below  the  dew  point  temperature  (e.g.  below 
60  °C),  not  only  much  more  sensible  heat  but  also  much  condens¬ 
ing  heat  will  be  recovered.  This  concept  is  called  as  exhaust  gas 
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Nomenclature 

do,  fli,  0-2 

polynomial  fitting  parameters 

Q 

thermal  power  (kW) 

A 

heat  exchanger  structure  parameter 

Re 

Reynolds  number 

C 

operation  cost 

T 

temperature  (°C) 

cp 

specific  heat  at  constant  pressure  (kj/kg  °C) 

TAGS 

a  kind  of  material  used  in  TEG 

CCHP 

combined  cooling  heating  and  power 

TEG 

thermoelectric  generator 

COP 

CHE 

coefficient  of  performance 
condensing  heat  exchanger 

zT 

thermoelectric  figure  of  merit  of  TEG 

CSR 

cost  saving  ratio 

Subscript 

d 

feature  size  of  exhaust  gas  pipe  (m) 

b 

boiler 

DAC 

double  effect  absorption  chiller 

c 

cold  side  of  TEG 

ESR 

energy  saving  ratio 

con 

conventional  separated  system 

F 

heat  exchange  area  (m2) 

dew 

dew  point  of  exhaust  gas 

f 

function  relationship 

el 

electricity 

h 

enthalpy  (kj/kg) 

exh 

exhaust  gas 

i 

index  number 

fuel 

fuel  consumption 

j 

index  number 

fw 

feed  water 

ICE 

internal  combustion  engine 

grid 

public  grid 

JHE 

jacket  water  heat  exchanger 

jw 

jacket  water 

k 

heat  transfer  coefficient  (kj/m2  °C) 

h 

hot  side  of  TEG 

L 

latent  heat  of  water  (kj/kg) 

loss 

heat  loss 

M 

molar  mass  (kg/mol) 

i 

inlet  flow 

m 

mass  flow  rate  (kg/s) 

m 

exhaust  gas  state  before  entering  CHE 

N 

number  of  TEG  modules 

0 

outlet  flow 

n 

molar  quantity 

rated 

rated  condition  of  ICE 

NG 

natural  gas 

sat 

saturation  state  of  water 

Nu 

Nusselt  number 

vap 

vapor  state  of  water 

P 

PEE 

electric  power  (kW);  pressure  (Pa) 
primary  energy  efficiency 

w 

water 

PLR 

part  load  ratio 

Greek  letter 

Pr 

Prandtl  number 

energy  price  ratio 

f 

Pri 

energy  price 

n 

efficiency 

deep-recovery  [25],  which  can  be  realized  by  adding  condensing 
heat  exchanger  (CHE).  In  fact,  the  condensing  heat  recovery  of 
exhaust  gas  has  been  widely  applied  in  condensing  boiler  in  Eur¬ 
ope,  and  the  efficiency  improvement  is  over  5%  [26-28].  If  deep- 
recovery  is  applied  in  CCHP  system,  the  primary  energy  efficiency 
will  be  obviously  enhanced. 

Basing  on  above  analyses,  a  framework  of  CCHP  system  is  pro¬ 
posed  for  power  generation,  refrigeration  and  domestic  hot  water 
production.  This  paper  proceeds  as  follow.  Firstly,  the  principle  of 
this  CCHP  system  is  described.  Secondly,  simulation  models  are 
developed  and  methodologies  are  introduced.  Thirdly,  system 
characteristics,  performance  are  discussed.  Finally,  some  main  con¬ 
clusions  are  summarized. 


2.  System  description 

Fig.  1  shows  the  schematic  diagram  of  this  CCHP  system.  Com¬ 
paring  with  some  commonly  used  CCHP  systems  [5— 7,9— 1 1 ,29— 
32],  DAC  is  applied  instead  of  single  effect  absorption  chiller  or 
adsorption  chiller,  TEG  is  applied  instead  of  exhaust  gas  heat 
exchanger,  and  CHE  is  added  to  further  recover  the  exhaust  gas 
waste  heat  by  deep-recovery  process.  Generally,  heating  load 
mainly  includes  house  heating  load  and  domestic  hot  water  load. 
In  this  system,  in  order  to  realize  exhaust  gas  deep-recovery,  the 
feed  water  temperature  should  be  low.  The  back  water  tempera¬ 
ture  from  house  heating  cycle  may  be  too  high  to  realize  exhaust 
gas  deep-recovery.  Therefore,  a  primary  assumption  is  that  only 
the  domestic  hot  water  load  is  satisfied  by  the  CCHP  system.  Also, 
it  can  be  understood  that  the  existing  hotels  and  residential 


buildings  in  the  north  of  China  have  been  equipped  with  central¬ 
heating  system  to  specially  meet  the  house  heating  load  demand. 

If  valve  1  is  turned  on  and  valve  2  is  turned  off,  the  system  will 
operate  in  combined  heating  and  power  mode  (CHP  mode).  In  CHP 
mode,  the  DAC  is  bypassed.  So,  there  is  no  cooling  output.  The 
exhaust  gas  from  ICE  successively  flows  through  TEG  and  CHE.  In 
TEG,  most  of  the  sensible  heat  of  exhaust  gas  is  absorbed  and  elec¬ 
tricity  is  generated.  In  CHE,  the  remained  small  part  of  sensible 
heat  and  some  condensing  heat  of  exhaust  gas  is  recovered  by 
deep-recovery  process.  Meanwhile,  the  feed  water  flows  through 
CHE,  TEG,  and  JHE  in  sequence  (1-2-4-5-6)  and  turns  into  hot 
water.  The  hot  water  is  entirely  consumed  as  domestic  hot  water 
by  user,  so  there  is  no  circulating  back  water.  Namely,  the  feed 
water  should  be  from  tap  water  which  is  cold.  In  CHE,  the  feed 
water  absorbs  the  low  temperature  waste  heat  by  deep-recovery 
process.  And  then  it  will  be  used  as  coolant  water  to  cool  the 
TEG.  In  TEG,  only  a  small  part  of  heat  can  be  converted  into  elec¬ 
tricity,  while  most  of  the  absorbed  heat  is  finally  taken  away  by 
the  coolant  water  in  the  cold  side  of  TEG.  So,  the  feed  water  can 
be  further  heated  by  the  cold  side  of  TEG.  After  flows  through 
the  TEG,  the  feed  water  flows  into  JHE  and  absorbs  the  whole  waste 
heat  in  jacket  water.  In  this  way,  the  cold  feed  water  is  heated  step 
by  step  and  finally  flows  into  the  buffer  tank. 

If  valve  1  is  turned  off  and  valve  2  is  turned  on,  the  system  will 
operate  in  combined  cooling,  heating  and  power  mode  (CCHP 
mode).  In  CCHP  mode,  the  TEG  is  bypassed.  So  there  is  no  electric 
output  from  TEG.  The  exhaust  gas  from  ICE  flows  into  DAC  for 
refrigeration,  and  then  it  flows  into  CHE  for  deep-recovery.  In 
DAC,  most  of  the  sensible  heat  of  exhaust  gas  is  recovered.  In 
CHE,  the  remained  small  part  of  low  temperature  sensible  heat 
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Fig.  1.  Schematic  diagram  of  CCHP  system  with  exhaust  gas  deep-recovery  and  thermoelectric  generator. 


and  some  condensing  heat  of  exhaust  gas  is  recovered  by  deep- 
recovery  process.  Meanwhile,  the  feed  water  flows  through  CHE 
and  JHE  in  sequence  (1 -3-5-6)  and  turns  into  hot  water  collected 
in  buffer  tank. 

3.  Modeling 

3.1.  ICE 


be  conducted  under  various  PLRICe .  The  rated  condition  of  DAC  is 
designed  based  on  the  exhaust  gas  condition  with  PLRICe  of  1. 

In  off-design  operation,  the  heat  transfer  performance  of  high 
pressure  generator  changes  with  the  mass  flow  rate  and  tempera¬ 
ture  of  exhaust  gas.  This  kind  of  relationship  can  be  simply  related 
by  [34]: 

kF  =  Am°xhTr  (3) 


The  performance  of  ICE  changes  with  the  part  load  ratio  of  ICE 
( PLRICe )•  The  PLRICe  represents  the  operating  state  which  is  defined 
by  the  electric  output  and  rated  electric  output  of  ICE,  as: 

PLRICe  =  Pei, iceI Pei, ice, rated-  Some  parameters,  such  as  electric  out¬ 
put,  waste  heat  of  jacket  water,  waste  heat  of  exhaust  gas,  exhaust 
gas  temperature,  fuel  consumption,  and  air  consumption  are  all 
related  with  PLRICE.  In  this  work,  the  relationships  between  these 
parameters  and  PLRICe  are  all  tested  by  experiment.  And  these 
experimental  parameters  constitute  the  start  point  of  the  following 
simulation. 

The  fuel  energy  consumed  by  ICE  is  finally  converted  into  elec¬ 
tricity,  waste  heat,  and  heat  loss.  So  the  energy  balance  in  ICE  is: 

Qfuel.lCE  =  Pel  T  Qjw  "T  Qexh  "T  Q./oss  (1) 


Tw.6  —  Tw 


>  '  ttl/w 


(2) 


The  hot  water  outlet  temperature  after  JHE  in  state  6  is  calcu¬ 
lated  as: where,  the  hot  water  temperature  in  state  5  can  be  got 
in  the  following  simulation. 


3.2.  DAC 


where  k  is  the  heat  transfer  coefficient,  A  is  a  parameter  related  to 
the  heat  exchanger  structure,  mexh  is  the  mass  flow  rate  of  exhaust 
gas,  r  is  a  factor  considering  the  influence  of  thermo-physical  prop¬ 
erties  of  exhaust  gas.  r  can  be  estimated  by: 

r  =  0.67s  -  0.27/  (4) 

where  s  and  l  are  the  powers  of  heat  conductivity  coefficient  and 
kinematic  viscosity  with  temperature,  respectively. 

The  total  mass  balance,  solute  mass  balance  and  energy  balance 
equations  for  each  component  of  DAC  can  be  respectively  repre¬ 
sented  as: 

Y/mi-J2mo  =  °  (5) 

y  m,Wi  -  m0w,-  =  0  (6) 

Q.+Y  m‘hi  -  Y  m°h°  = 0  (?) 

where  the  subscript  i  and  o  mean  inlet  flow  and  outlet  flow, 
respectively. 

The  heat  transmission  in  each  component  of  DAC  can  be  calcu¬ 
lated  by: 


The  DAC  applied  in  this  work  is  a  lithium  bromide  water  (LiBr- 
H20)  type  with  the  high  pressure  generator  and  low  pressure  gen¬ 
erator  in  parallel.  The  condensing  heat  of  refrigerant  vapor  from 
high  pressure  generator  is  used  to  power  the  low  pressure  genera¬ 
tor  firstly,  and  then  part  of  its  sensible  heat  is  further  recovered  by 
low  pressure  heat  exchanger.  This  type  of  structure  can  be  referred 
to  [33].  In  simulation,  the  high  pressure  generator  of  DAC  is  pow¬ 
ered  by  exhaust  gas  directly,  and  the  DAC  is  cooled  by  a  cooling 
tower.  As  powered  by  the  waste  heat  of  exhaust  gas  from  ICE, 
the  performance  of  DAC  depends  on  PLRICe .  So,  the  simulation  will 


Q  =  kFAT  (8) 

where  AT  is  the  logarithmic  mean  temperature  difference.  In  high 
pressure  generator,  kF  can  be  calculated  by  Eq.  (3).  In  any  other 
component,  the  temperature  fluctuation  is  small  and  kF  can  be 
assumed  as  constant  which  can  be  calculated  by  simulation  in  rated 
condition. 

All  the  thermodynamic  properties  of  LiBr-H20  solution  used  in 
simulation  come  from  results  of  [35].  Some  other  commonly 
applied  equations  and  assumptions  in  simulation  can  be  referred 
to  [33]  and  [36].  By  simulation,  the  exhaust  gas  outlet  temperature 
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Texh,m,  refrigeration  COP,  and  some  other  results  under  various 
PLRICe  can  be  calculated. 

The  heat  consumption  by  DAC  is  calculated  by: 

Q.DAC  =  Cp.exhfflexh  (Texh,  i  ~  Texh,m)  (9) 

The  cooling  output  of  DAC  is  calculated  by: 

Qc  —  COPqac  •  Qdac  (10) 


3.3.  TEG 


TEG  is  a  device  that  converts  thermal  power  into  electric  power 
based  on  the  principle  of  Seebeck  effect  [37].  If  there  is  a  temper¬ 
ature  difference  between  the  hot  side  and  cold  side  of  TEG,  there 
will  be  electric  current  through  the  TEG.  TEG  usually  consists  of 
many  TEG  modules.  The  transport  characteristics  of  a  TEG  module 
can  be  represented  by  the  thermoelectric  figure  of  merit  (zT).  For  a 
given  material,  the  efficiency  of  a  TEG  module  depends  of  the  tem¬ 
perature  condition,  which  can  be  calculated  as  38]: 


AT 

t,~Th  +  273.15 


Vl  +  zT  -  1 


\/T+zT  + 


rc+273.15 
Th+ 273.15 


(11) 


where  zT  is  the  thermoelectric  figure  of  merit  calculated  at  mean 
temperature  between  hot  side  and  cold  side. 

Similar  to  a  heat  engine,  the  theoretical  maximum  electric  effi¬ 
ciency  of  a  TEG  module  is  limited  by  the  Carnot  efficiency,  namely 
ATl(Th  +  273.15).  A  larger  zT  will  help  to  increase  the  efficiency 
under  a  given  temperature  condition.  In  this  work,  TAGS,  which 
has  been  successfully  used  in  long-life  thermoelectric  generators, 
is  used  as  the  TEG  module  material.  Snyder  et  al.  presented  a  func¬ 
tion  curve  for  zT  of  TAGS  with  temperature  38].  Herein,  this  func¬ 
tional  relationship  is  approximately  described  by  a  polynomial 
fitting  equation  as: 

zT(T)  =  a2  ■  T2  +  cii  T  +  a0  (12) 


where  a2,  cq  and  a0  are  -6.879e-6,  6.443e-3,  and  -0.3498, 
respectively. 

Technically,  the  TEG  modules  usually  distribute  around  the 
exhaust  gas  pipe  as  some  products  show  [22].  The  hot  side  of  a 
TEG  module  contacts  with  the  exhaust  gas.  The  Nusselt  number 
of  heat  transfer  can  be  estimated  by  the  Dittus-Boelter  formula 

as  [39]: 

Nu  =  0.023ReosPr°3  (13) 


As  the  hot  side  temperature  changes  along  the  exhaust  gas  flow, 
the  performance  of  a  TEG  module  depends  on  the  local  condition. 
The  heat  absorbed  by  the  hot  side  of  a  TEG  module  is: 

Q TEG,i  =  ^F(TeXh,<  ~  ThJ)  (14) 

where  i  represents  ith  TEG  module,  i  =  1,. . .,  N.  The  total  electric  out¬ 
put  of  TEG  should  be  calculated  as: 


Pteg  = 


(15) 


In  turn,  the  efficiency  of  TEG  device  should  be  calculated  as: 

Pteg 


Vteg  — 


Cp,exhWlexh{Texh,i  P  exh,m ) 


(16) 


Then,  the  heat  absorbed  by  the  coolant  water  in  the  cold  side  of 
TEG  can  be  calculated  as: 


Qteg,c  —  4c(Qteg,i i  -  Pteg)  (17) 

where  heat  recovery  efficiency  rjc  in  cold  side  of  TEG  is  assumed  as 
0.9  considering  heat  loss  in  TEG. 


The  feed  water  absorbs  the  heat  from  the  cold  side  of  TEG  and 
then  the  water  temperature  increases  from  state  2  to  state  4: 


Tw,  4  —  TW2  + 


Q.TEG,c 

Cp,w  ■  aifw 


(18) 


3.4.  Exhaust  gas  deep-recovery 

For  a  given  combustion  condition,  including  fuel  composition 
and  air  consumption  rate,  the  exhaust  gas  composition  can  be  cal¬ 
culated.  Once  the  exhaust  gas  composition  is  known,  the  conden¬ 
sation  characteristics  can  be  analyzed.  Herein,  it  is  assumed  that 
the  fuel  is  completely  burned.  And  the  exhaust  gas  is  assumed  as 
ideal  gas,  which  consists  of  02,  C02,  N2  and  H20. 

The  dew  point  temperature  Tdew  of  exhaust  gas  depends  on  the 
partial  pressure  of  gaseous  water  Pvap. 

For  ideal  gas,  Pvap  can  be  calculated  by  the  mole  fraction  as: 

Pvap  —  Pexh  '  H-vap/tlexh  (19) 

where  nvap  is  the  molar  quantity  of  gaseous  water,  nexh  is  the  molar 
quantity  of  exhaust  gas,  and  both  of  them  can  be  calculated  accord¬ 
ing  to  the  chemical  equation  of  combustion.  The  higher  the  gaseous 
water  content  is,  the  higher  the  dew  point  temperature  will  be.  The 
relationship  between  dew  point  temperature  and  partial  pressure  of 
gaseous  water  can  be  associated  by: 

T dew  —  fvap  (P vap  )  (20) 

where  fvap  represents  the  thermodynamic  function  for  gaseous 
water,  which  can  be  calculated  by  REFPROP  [40]. 

The  gaseous  water  content  before  CHE  can  be  calculated  by: 

ttlvapyn  —  Hyap  '  M vap  (21) 

where  Mvap  is  the  molar  mass  of  water. 

And  the  gaseous  water  content  after  CHE  can  be  calculated  by: 

Hi vap, o  —  “  “  '  fleXh,o  '  M Vap  (22) 

Pexh,o 

where  Psat,o  is  the  saturation  pressure  of  gaseous  water  correspond¬ 
ing  to  the  dew  point  of  exhaust  gas  after  CHE.  Then  the  mass  of  con¬ 
densed  water  mdew  can  be  got  by: 

dew  —  iilyap.m  ~  HI  vap, o  (23) 

The  recoverable  waste  heat  in  the  process  of  deep-recovery  is 
calculated  by: 

4 

0.CHE  =  '  (fy>  “  ho)  +  mdew  '  L  (24) 

J=1 

where  j  =  1,  2,  3,  4  corresponds  to  02,  C02,  N2,  and  H20  in  gaseous 
state,  respectively.  L  is  the  latent  heat  of  phase  change  of  water. 

The  feed  water  is  firstly  heated  by  CHE  and  then  the  water  tem¬ 
perature  increases  to  state  2  (or  3): 

Tw, 2  —  Tw>i  +  Qche/ {Cp,w  ■  aifw)  (25) 


3.5.  Performance  evaluation 


Primary  energy  efficiency  (PEE)  is  used  to  describe  the  total  effi¬ 
ciency  of  CCHP  system  based  on  the  first  law  of  thermodynamics.  It 
represents  the  ratio  of  the  whole  energy  output  to  primary  energy 
consumption  (lower  heating  value)  [5]. 


p££  _  Pei  +  Qc  +  Qh 

Qf uel, CCHP 


(26) 


Primary  energy  saving  ratio  (ESP)  and  cost  saving  ratio  ( CSR )  are 
widely  used  to  evaluate  the  performance  of  CCHP  system  compared 
with  conventional  separated  system  [41].  The  ESR  is  defined  as: 
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ESR  = 


Qfuel  yon  Qfuel, CCHP 

Qfuel  yon 


(27) 


Typically,  conventional  separated  system  is  assumed  as:  the 
electric  load  is  satisfied  by  the  public  grid,  the  cooling  load  is  sat¬ 
isfied  by  an  electric  chiller,  and  the  heating  load  is  satisfied  by  a  gas 
boiler.  So,  the  fuel  consumption  of  conventional  separated  system 
is  calculated  by: 


Qfuelyon 


Pel  +  Qc/COPel 

yhiller 


CSR  = 


Vgrid  4  b 

Similarly,  the  CSR  is  defined  as: 

CCon  ~  CcCHP 


(28) 


(29) 


The  operation  cost  is  defined  according  to  the  cost  by  energy 
consumption: 

C  —  Pel.gridPriel  +  {Qfuel, b  +  Qfuel, i ICe)  P^fuel  (30) 


where  Qfuei,b  and  Pei,grid  are  both  0  for  CCHP  system,  Q/uel/C£  is  0  for 
conventional  separated  system. 

So,  Eq.  (29)  can  be  calculated  as: 


prn  _  (^el-grid  +  £  '  Qfuel, b)con  (j  '  Qfuel, ICe)CchP 

(Pel, grid  +  i  '  Qfuel, b) con 


(31) 


where  £  is  the  energy  price  ratio  defined  as  £  =  Prifuei/Priei  [40]. 


4.  Results  and  discussion 

4.1.  Exhaust  gas  deep-recovery 

A  16  kW  ICE  is  tested  and  some  related  experimental  results  are 
listed  in  Table  1.  The  ICE  is  tested  from  idling  condition  to  full  load 
condition,  with  PLRICe  changes  from  0  to  1.  The  exhaust  gas  tem¬ 
perature  increases  with  PLRICe,  and  it  reaches  670.4  °C  at  full  load 
condition. 

Table  2  shows  the  composition  of  NG.  The  content  of  CH4  is 
97.92%.  So,  there  will  be  much  gaseous  water  contained  in  the 
exhaust  gas  of  NG.  In  deep-recovery,  most  of  the  gaseous  water 
can  be  condensed  and  the  condensing  heat  can  be  recovered  by 
cold  feed  water.  Basing  on  the  NG  consumption  rate,  air  flow  rate, 
and  fuel  composition,  the  dew  point  and  condensing  heat  of 
exhaust  gas  can  be  calculated.  As  shown  in  Fig.  2,  the  condensing 
heat  changes  from  1.37  to  4.4  kW  when  the  PLRICe  changes  from 
0.1  to  1.  The  corresponding  dew  point  temperature  of  exhaust 
gas  changes  from  54.7  to  57.1  °C,  which  is  a  little  lower  than  the 
result  calculated  by  Che  et  al.  [26].  As  the  air  factor  in  their  studied 
condition  is  lower  than  this  work,  they  concluded  that  the  dew 
point  temperature  is  generally  56-60  °C. 

The  recoverable  waste  heat  in  deep-recovery  includes  sensible 
heat  and  condensing  heat.  Fig.  3  shows  the  recoverable  waste  heat 


Table  1 

Experimental  results  of  a  1 6  kW  ICEa. 


PLR,ce 

NG  consumption  rate 

(N  m3/h) 

Air  flow  rate 

(N  m3/h) 

Exhaust  gas 
temperature  (°C) 

0 

2.01 

26.44 

439.46 

0.1 

2.32 

31.60 

473.9 

0.2 

2.56 

34.77 

480.95 

0.3 

2.88 

39.29 

512.36 

0.4 

3.25 

43.28 

523.39 

0.5 

3.75 

48.12 

548.93 

0.6 

4.25 

52.11 

567.5 

0.7 

4.78 

57.78 

594.87 

0.8 

5.4 

64.29 

603.3 

0.9 

6.07 

73.95 

627.95 

1 

6.76 

84.67 

670.4 

a  The  type  of  ICE  is  Cummins  GGDB16. 


of  ICE  and  total  recoverable  waste  heat  increasing  rate  by  deep- 
recovery.  As  the  PLRICe  increases  from  0.1  to  1,  the  waste  heat  of 
jacket  water  varies  from  11.2  kW  to  22.3  kW,  and  the  waste  heat 
of  exhaust  gas  (>100  °C)  varies  from  4.97  to  14.2  kW.  If  the  exhaust 
gas  is  further  cooled  to  40  °C  by  deep-recovery  (It  has  been  realized 
in  a  preliminary  test  by  20  °C  feed  water),  much  more  waste  heat 
can  be  recovered  as  the  white  columns  show.  It  makes  the  recover¬ 
able  exhaust  gas  waste  heat  and  total  recoverable  waste  heat 
increase  by  41-46%  and  13-16%,  respectively.  So,  exhaust  gas 
deep-recovery  can  greatly  improve  the  energy  utilization  efficiency. 

4.2.  Performance  in  CCHP  mode 

Fig.  4  shows  the  part  load  ratio  performance  of  system  in  CCHP 
mode.  As  the  TEG  does  not  operate,  all  the  electric  output  comes 
from  ICE  and  the  maximum  electric  output  is  16  kW.  The  cooling 
output  increases  from  6.14  kW  to  17.56  kW  while  the  heating  out¬ 
put  increases  from  12.9  kW  to  30.31  kW.  It  means  that,  even  in 
CCHP  mode,  the  heating  output  is  obviously  much  more  than  the 
cooling  output.  So,  this  kind  of  CCHP  system  is  suitable  for  users 
with  large  domestic  hot  water  demand.  Analysis  shows  that, 
13.3-19.5%  fuel  energy  (lower  heating  value)  is  lost  to  the  environ¬ 
ment  in  this  system,  but  deep-recovery  of  exhaust  gas  contributes 
some  other  condensing  heat  which  does  not  count  in  the  lower 
heating  value.  In  addition,  the  COP  of  DAC  is  great  than  1,  which 
also  magnifies  the  thermal  output.  So,  the  average  PEE  reaches 
0.944  as  shown  in  Fig.  4.  The  theoretical  PEE  can  be  expected  to 
be  over  1  if  condensing  heat  is  recovered  as  much  as  possible. 

Fig.  5  shows  the  influence  of  PLRICe  on  performance  of  DAC.  Both 
the  exhaust  gas  outlet  temperature  after  DAC  and  thermal  COP  are 
investigated.  In  the  rated  condition  with  PLRICe  of  1,  the  exhaust  gas 
outlet  temperature  is  designed  at  170°C.  If  PLRICe  decreases  to  0, 
namely  in  idling  condition,  the  exhaust  gas  outlet  temperature  will 
decreases  to  107  °C.  Simultaneously,  the  COP  of  DAC  presents  a 
trend  of  parabolic  curve.  The  COP  of  DAC  reaches  maximum  when 
the  PLRICe  is  less  than  1,  which  means  that  the  best  heat  and  mass 
transfer  performance  of  DAC  is  not  in  full  load  condition.  Zheng 
et  al.  [34]  and  Wu  et  al.  [42]  have  also  respectively  pointed  out  that 
the  best  performance  of  DAC  and  ICE  may  not  in  full  load  condition. 

The  hot  water  to  the  buffer  tank  represents  the  whole  heating 
output.  Obviously,  the  hot  water  outlet  temperature  is  important 
for  user.  As  shown  in  Fig.  6,  it  is  the  contour  map  of  hot  water  out¬ 
let  temperature  in  CCHP  mode,  considering  the  influence  of  PLRICe 
and  feed  water  flow  rate.  The  orange  zone,  which  is  named  as  fea¬ 
sible  operating  zone,  represents  hot  water  outlet  temperature 
ranging  between  45  °C  and  85  °C.  These  two  temperature  values 
are  assumed  artificially.  A  lower  temperature  may  be  unacceptable 
by  user.  A  higher  temperature  may  influence  the  cooling  of  jacket 
water  cycle  and  damage  the  ICE.  For  a  given  PLR1Ce,  which  means 
the  waste  heat  is  given,  the  hot  water  outlet  temperature  increases 
as  the  feed  water  flow  rate  decreases.  For  a  given  feed  water  flow 
rate,  the  hot  water  outlet  temperature  increases  with  PLRICe •  It  can 
be  found  that  the  hot  water  outlet  temperature  will  be  always 
within  the  assumed  temperature  range  if  the  feed  water  flow  rate 
is  set  between  0.4  m3/h  and  0.44  m3/h. 

4.3.  Performance  in  CHP  mode 

Fig.  7  shows  the  part  load  ratio  performance  of  system  in  CHP 
mode.  As  the  DAC  does  not  operate,  there  is  no  cooling  output. 
The  electric  output  comes  from  the  ICE  and  TEG  simultaneously. 
Both  the  electric  output  and  heating  output  increase  with  the 
PLRICe •  The  electric  output  of  system  reaches  17.18  kW  when  PLRICe 
is  1,  in  which  16  kW  comes  from  ICE  and  the  extra  1.18  kW  comes 
from  TEG.  The  heating  output  comes  from  the  released  heat  of  the 
cold  side  of  TEG,  the  waste  heat  of  exhaust  gas  in  deep-recovery, 
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Table  2 

Composition  of  NG. 


Contents 

co2 

n2 

ch4 

c2h6 

c3h8 

C4H10 

c5h12 

c6 

Percentage  (%) 

0.657 

0.72 

97.92 

0.561 

0.058 

0.037 

0.013 

0.034 

3 

crq 

tr 


5 


1.40 

1.38 

1.36 

1.34 

1.32 

1.30 

1.28 

1.26 

1.24 

1.22 

1.20 


Part  load  ratio  of  ICE 

Fig.  5.  Influence  of  PLRICe  on  performance  of  DAC. 
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Fig.  3.  Recoverable  waste  heat  of  ICE  and  total  recoverable  waste  heat  increasing 
rate  by  deep-recovery. 
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Fig.  4.  Part  load  ratio  performance  of  system  in  CCHP  mode. 
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Fig.  6.  Contour  map  of  hot  water  outlet  temperature  in  CCHP  mode. 


and  the  waste  heat  of  jacket  water.  It  can  be  found  that  the  average 
PEE  reaches  0.908.  This  value  agrees  with  the  result  reported  by  Fu 
et  al.  [43].  In  their  experiment,  by  recovering  the  jacket  water 
waste  heat,  sensible  heat  of  exhaust  gas  and  condensing  heat  of 
exhaust  gas  for  heating,  the  PEE  reaches  about  90-94%. 

Fig.  8  shows  the  influences  of  PLRICe  on  electric  output  of  TEG 
and  electric  efficiency  of  CCHP  system.  It  can  be  found  that  the 
electric  output  of  TEG  increases  with  PLRICe .  The  minimum  electric 
output  of  TEG  is  0.23  kW,  which  occurs  when  the  ICE  operates  in 
idling  condition  with  PLRICe  of  O.The  maximum  electric  output  of 
TEG  is  1.18  kW,  which  occurs  when  the  ICE  operates  in  full  load 
condition  with  PLRICe  of  1.  As  a  demonstration,  Bass  et  al.  have  real¬ 
ized  1  kW  power  output  by  installing  a  TEG  device  outside  the 
exhaust  gas  pipe  of  a  truck  [21  ].  As  TEG  is  applied,  the  electric  effi¬ 
ciency  of  CCHP  system  is  improved.  When  the  PLRICe  is  1,  the  elec¬ 
tric  efficiency  is  improved  from  0.237  to  0.254  with  increasing  rate 
of  7.4%. 

Fig.  9  shows  the  influence  of  PLRICe  on  performance  of  TEG.  Both 
the  exhaust  gas  outlet  temperature  after  TEG  and  electric  efficiency 
of  TEG  are  investigated.  In  the  rated  condition  with  PLRICe  of  1 ,  the 
exhaust  gas  outlet  temperature  is  designed  to  be  170  °C.  Both  the 
exhaust  gas  outlet  temperature  and  electric  efficiency  increases 
with  PLRICe •  If  PERice  decreases  to  0,  namely  the  ICE  operates  in 
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Part  load  ratio  off  ICE 
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Fig.  7.  Part  load  ratio  performance  of  system  in  CHP  mode. 


Fig.  9.  Influence  of  PLRICe  on  performance  of  TEG. 
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Fig.  8.  Influences  of  PLRICe  on  electric  output  of  TEG  and  electric  efficiency  of  CCHP 
system. 

idling  condition,  the  outlet  temperature  will  decrease  to  91  °C.  The 
electric  efficiency  of  TEG  can  reach  0.087  and  0.044  with  PLRICe  of  1 
and  0,  respectively.  In  comparison,  the  efficiencies  are  reported  at 
0.09  and  0.065  based  on  module  material  of  EC021  for  the  same 
temperature  differences,  respectively  [44].  So,  the  performance  of 
TEG  is  greatly  influenced  by  PLRICe .  A  higher  PLRICe  means  higher 
exhaust  gas  temperature  and  flow  rate  from  ICE,  which  can 
improve  the  performance  of  TEG. 

Fig.  10  shows  the  influence  of  feed  water  flow  rate  on  the  per¬ 
formance  of  TEG  with  PLRice  of  1 .  A  larger  feed  water  flow  rate  can 
help  to  cool  the  cold  side  of  TEG  and  improve  the  performance  of 
TEG.  So,  the  electric  efficiency  of  TEG  increases  with  the  flow  rate. 
Simultaneously,  the  water  outlet  temperature  of  TEG  decreases  as 
the  feed  water  flow  rate  increases.  So,  a  larger  flow  rate  will  result 
in  a  lower  hot  water  outlet  temperature.  There  should  be  a  feasible 
range  of  feed  water  flow  rate. 

As  Figs.  1 1  and  12  show,  they  are  the  contour  maps  of  hot  water 
outlet  temperature  in  the  transitional  season  condition  and  winter 
season  condition,  respectively.  According  to  the  typical  tap  water 
temperature,  the  feed  water  temperatures  in  the  transitional  sea¬ 
son  and  winter  season  are  assumed  at  12  °C  and  5  °C,  respectively. 
As  the  orange  feasible  operating  zone  shows,  too  large  or  too  small 
feed  water  flow  rate  will  shrink  the  feasible  operating  range  of 
PLRICe •  For  the  transitional  season  condition,  a  flow  rate  of 
0.45-0.5  m3/h  will  be  the  most  adaptive.  For  the  winter  season 
condition,  a  flow  rate  of  0.38-0.47  m3/h  can  be  considered.  If  the 
hot  water  should  be  temporarily  stored  in  the  buffer  tank,  the  flow 
rate  can  be  controlled  as  low  as  possible  to  increase  the  stored  hot 
water  temperature. 


Fig.  10.  Influence  of  feed  water  flow  rate  on  the  performance  of  TEG  (PLRICe=  1). 

4.4.  System  evaluation 

Fig.  13  shows  the  influence  of  PLRICe  on  ESR.  As  the  black  curves 
show,  system  operations  without  TEG  and  CHE  are  used  for  com¬ 
parisons.  The  ESRs  of  CCHP  mode  and  CHP  mode  are  similar  to  each 
other.  But,  there  is  obvious  difference  for  ESR  between  systems 
with  and  without  TEG  and  CHE.  Especially  in  CHP  mode  with  PLRICe 
of  1,  the  ESR  is  improved  from  0.216  to  0.304  by  applying  TEG  and 
CHE.  When  the  PLRICe  is  over  0.5,  the  ESR  of  CCHP  system  is  about 
0.28  while  it  is  only  about  0.21  without  TEG  and  CHE.  So,  as  TEG 
and  CHE  are  applied,  much  more  energy  can  be  saved.  Once  the 
PLRICe  is  below  0.5,  the  system  performances  dramatically  get 
worse.  This  similar  trend  also  occurs  to  CCHP  system  without 
TEG  and  CHE.  So,  CCHP  systems  are  recommended  to  be  operated 
with  high  PLRICe • 

Fig.  14  shows  the  influence  of  PLRICe  on  CSR.  In  CCHP  mode,  the 
CSR  of  system  is  about  0.04-0.09  higher  than  system  without  TEG 
and  CHE.  In  CHP  mode,  it  is  about  0.05-0.12.  For  CHP  mode  with 
PLRICe  of  1,  the  CSR  of  system  reaches  0.417  while  it  is  0.349  with¬ 
out  TEG  and  CHE.  So,  by  applying  TEG  and  CHE,  the  system  opera¬ 
tion  can  become  much  more  economical. 

Further  considering  a  floating  energy  price,  Fig.  15  shows  the 
sensitivity  analysis  of  energy  price  ratio  on  CSR.  The  theoretical 
maximum  CSR  is  1  with  energy  price  ratio  of  0,  namely  the  fuel 
price  is  0  or  the  electric  price  is  infinitely  great.  A  higher  fuel  price 
means  a  worse  economical  performance  of  CCHP  system  operation. 
So,  the  CSR  decreases  as  the  energy  price  ratio  increases.  The  con¬ 
tinuous  curves  represent  the  performances  of  CCHP  system,  and 
the  dash  curves  represent  the  performances  of  CCHP  system  with¬ 
out  TEG  and  CHE.  The  red  dash  lines  indicate  the  present  local  case 
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Fig.  11.  Contour  map  of  hot  water  outlet  temperature  in  the  transitional  season 
condition. 
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Fig.  12.  Contour  map  of  hot  water  outlet  temperature  in  winter  season  condition. 
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Fig.  13.  Influence  of  PLRICe  on  ESR. 

with  energy  price  ratio  of  0.25.  The  red  points  marked  on  the  X  axis 
show  the  energy  price  ratios  with  CSR  of  0.  They  are  named  as  crit¬ 
ical  CSRs.  If  TEG  and  CHE  are  not  applied,  the  critical  CSRs  in  CCHP 
mode  and  CHP  mode  are  0.49  and  0.58,  respectively.  If  TEG  and 
CHE  are  both  applied,  the  critical  CSRs  in  CCHP  mode  and  CHP 
mode  get  much  higher  and  they  are  0.62  and  0.84,  respectively. 
Obviously,  the  economical  adaptability  is  greatly  improved. 

As  TEG  and  CHE  are  incorporated,  the  initial  investment  should 
be  higher.  Table  3  shows  the  investments  of  related  devices  and 
system  integration  fees.  Basing  on  Table  3,  Fig.  16  shows  the 


Part  load  ratio  of  ICE 


Fig.  14.  Influence  of  PLRICe  on  CSR. 


Table  3 

Prices  and  investments  of  every  involved  device. 


Items 

Price  ($/ 
kW) 

Capacity 

(kW) 

Investment 

(S) 

DAC 

900 

18 

16,200 

ICE-G 

800 

16 

12,800 

TEG 

2000 

1.2 

2400 

CHE 

70  [45] 

9 

630 

JHE 

10 

23 

230 

System  integration  without  TEG 
and  CHE 

- 

- 

7000a 

System  integration  with  TEG  and 
CHE 

" 

" 

8000b 

a  Considering  the  cost  of  control  system,  pipes,  valves,  and  so  on. 
b  Including  the  increased  investment  of  exhaust  gas  heat  exchanger  needed  by 
TEG. 


CCHP  system  with  TEG  and  CHE 
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Fig.  16.  Investments  comparison  between  systems  with  and  without  TEG  and  CHE. 
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investments  comparison  between  systems  with  and  without  TEG 
and  CHE.  It  can  be  found  that,  the  total  investment  is  mainly 
caused  by  the  ICE  and  DAC.  Although  the  present  TEG  price  for 
per  kW  is  very  high,  the  required  capacity  is  small.  So,  the  total 
investment  increment  of  system  is  about  11.1%. 

5.  Conclusions 

A  CCHP  system  is  proposed  for  power  generation,  refrigeration, 
and  domestic  hot  water  production.  The  operating  characteristics 
under  various  working  conditions  are  simulated  and  discussed. 
This  kind  of  CCHP  system  shows  an  obviously  better  performance 
with  some  increment  of  total  investment.  Some  main  conclusions 
are  summarized  as  follows. 

Firstly,  as  deep-recovery  of  exhaust  gas  is  applied,  the  total 
recoverable  waste  heat  of  CCHP  system  is  improved  by  13-16%. 
The  PEE  of  CCHP  system  reaches  0.944  and  0.908  in  CCHP  mode 
and  CHP  mode,  respectively. 

Secondly,  as  TEG  is  applied  to  generate  electricity  by  waste  heat, 
the  CCHP  system  can  realize  a  higher  electric  efficiency.  When 
PLRICe  is  1,  the  electric  efficiency  is  improved  from  0.237  to  0.254 
with  increasing  rate  of  7.4%. 

Thirdly,  there  is  a  feasible  operating  zone  in  order  to  control  the 
hot  water  outlet  temperature  between  45  °C  and  85°C.  According 
to  the  contour  maps  of  hot  water  outlet  temperature,  the  feasible 
controlling  strategy  of  ICE  and  feed  water  flow  rate  can  be  found. 

Fourthly,  the  system  performance  improvement  is  obvious  by 
applying  TEG  and  CHE.  For  CHP  mode  with  PLRICE  of  1,  the  ESR  is 
improved  from  0.216  to  0.304,  and  the  CSR  is  improved  from 
0.349  to  0.417.  In  addition,  the  economical  adaptability  is  also 
greatly  improved,  and  the  total  investment  increment  of  CCHP  sys¬ 
tem  is  11.1%. 
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